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ABSTRACT 

We present a study of the kinematic properties of the ionized gas in the dom- 
inant giant Hll region of the well known Hn galaxy: II Zw 40. High spatial 
and spectral resolution spectroscopy has been obtained using IFU mode on the 
GMOS instrument at Gemini-North§ telescope. The observations allow us to 
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obtain the Ha intensity map, the radial velocity and velocity dispersion maps 
as well as estimate some physical conditions in the inner region of the starburst, 
such as oxygen abundance (O/H) and electron density. We have used a set of 
kinematics diagnostic diagrams, such as the intensity versus velocity dispersion 
(7-cr), intensity versus radial velocity (I-V) and V-a, for global and individual 
analysis in sub-regions of the nebula. We aim to separate the main line broad- 
ening mechanisms responsible for producing a smooth supersonic integrated line 
profile for the giant Hn region. Bubbles and shells driven by stellar winds and 
possibly supernovae, covering a large fraction on the face of the nebula, are identi- 
fied on scales larger than 50 pc. We found that unperturbed or "free from shells" 
regions showing the lowest o values (~ 20 km s _1 ) should be good indicators 
for the Ograv component in II Zw 40. The brightest central region (R ~ 50 pc) 
is responsible for a derived from a single fit to the integrated line profile. The 
dominant action of gravity, and possibly unresolved winds of young (<10 Myr) 
massive stars, in this small region should be responsible for the characteristic 
Ha velocity profile of the starburst region as a whole (er = 32-35 kms -1 ). Our 
observations show that the complex structure of the interstellar medium of this 
galactic scale star-forming region is very similar to that of nearby extragalactic 
giant Hn regions in the Local Group galaxies. 

Subject headings: Galaxies: individual II Zw 40, Galaxy: kinematics and dynam- 
ics 



1. Introduction 



Hll galaxies (HllGs) are dwarf galaxies undergoing a burst of star formation. Their 
extensive star-forming regions, composed of an ensemble of Super Star Clusters (SSC) and 
giant Hll regions (GHnRs), dominate the optical emission. HllGs represent the simplest 
examples of the starburst phenomenon occurring on galactic scales. II Zw 40 is one of the 
most famous HnG showing an intense starburst in its central region. Its optical spectrum 
is dominated by very intense emission lines of H, He, [O n], [O ill], [N n], [S n] and [Ne 
ill] superposed on a faint blue continuum, indicati ng that the present rate of star formation 



(SF) is much higher than the historical average (ISargent feSearlgJ Il970l ; ISearle &: Sargent 



19721 : iTerlevich et al.lll99ll : IWalsh fc Rovill993l : IKehrig et al 



200 



Eh 



Zw 40 has a subsolar 



oxygen abundance of 12+log(O/H)=8.07 from the data of IKehrig et al.l (120041 ) . typical of 
this class. The star formation (SF) history of this galaxy has been the object of several 
studies, ranging from the hypothesis that the star-forming episode may be producing a first 
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generation of stars, to a scenario in which the history of SF consists of a conti nuous series 



of sho rt bursts over its lifetime. Near infrared spectra have been obtained by IVanzi et al. 



( 119961 ). who showed that II Zw 40 has low supernova rate and that no nuclear starburs t, as 
powerful as the presen t one, has occurred in the past Giga-year. iKunth fc Sargentl (119811 ) and 
Vacca fc Contil (119921 ) found Wolf-Rayet (WR) fea tures in the optica l spectrum, confirming 
that the central region is very young (3-6 Myr). iBeck et al.l (120021 ) found that the bright 
central region concentrate 75% of the thermal free-free emission and that the 75 pc r egion 
is formed by four supernebula each powered by ~ 600 O stars. lUlvestad et al.l (120071 ) also 
found no radio supernovae with powers greater than Cas A. 

The initial trigger of the present SF episode is also a question that remains unanswered 



for this class of dwarf starbursts since most are isolated systems (ITelles fc Terlevichl 11995 



Telles fc Maddoxll2000l ). If not triggered by external agents (i.e. inter actions), star formatio n 
is caused by internal processes manifesting in a stochastic manner (jPelupessy et al.l 12004 ). 
This may be the c ase for the subsample of low l umin osity, compact objects, classified as type 
II HII galaxies by lTelles. Melnick fc Terlevichl (119971 ). with no signs of morphological distur- 
bances. However, the more luminous HII galaxies, classified as type I, do show some signs of 
extensions, fuzz, or tails in their outer envelopes, suggesting tidal origin. Over the years, only 
a few studies have been carried out to determine the overall dynamics of II Zw 40. Hi veloc- 



ity mapping reveals complex structure with a subtle velocity gradient (IBrinks fc Kleinlll988 



van Zee et al.lll998l ). The southeast and northwest Hi tails have reversal rotation, suggesting 
that most of the material is falling back toward the dynamical center. The best scenario 
for the formation of II Zw 40 seems to be the result a collision between two gas-rich dwarf 



galax ies, and the current starburst occurs in the overlapping merging region ([Baldwin et al. 
1982h . 



In the early 70's, ISmith fc Weedmanl (ll970al 1972) first discovered broad [OIII] and 
Balmer lines in GHllRs. These lines have widths which are broader (i.e. line of sight velocity 
dispersion a ~ 13 — 25 kms" 1 ) than those found in normal Hll regions in the Galaxy, e.g. a 
~ 5.4 km s _1 in Orion (jMunchlll958l ; ISmith &: Weedmanlll970bl ). These high velocity widths 
imply gas motions a few times the speed of sound. The early detailed studies suggested 
that winds from WR star s should contribute to the total kinetic energy in mass motions 
(ISmith &: Weedmanl 119721 ). and that su pernova remnants (SNR s) could play a significant 
role in the energetic balance of GHllRs (jSkillman fc Balicklll984j ). 



In the meantime, iMelnickl (Il977l) showed a correlati o n bet ween the diameter and the 
global Ha profile widths of GHllRs. iTerlevich k Melnickl Jl98lh found L(H/3) oc a~ 4 which 
led them to propose a gravitational model for the origin of the supersonic motions, due to the 
close similarities with the parametric relations found for systems in dynamical equilibrium, 



-4- 



such as elliptical galaxies, g lobular cluster s and bulges of spiral galaxies. The first extrapo- 
lation to HnGs was done by lMelnick et al.l (119881 ). where a significant sample of galaxies was 
observed, including II Zw 40. Recent work also confirms the existence of the L-a relation for 
HnGs, suggest ing that GHllRs and the star-forming regions in HllGs belong to a single fam- 



ily of objects (IRozas et al.l 120061 ; iBosch et al.l 12002 ; iTelles et al.ll200ll ; iFuentes-Masip et al. 
2000h . 



Several line broadening mechanisms have been pro posed t o inter pret the supersonic 
widths measured in the spectra of GHllRs and HnGs. iDysonl (119791 ) proposes that line 
profiles consist of two components: one broad line component of the shocked region, and a 
narrower component originated in the neighboring region which is ionized by UV photons 
escaping from the broad line region. The broad line emission is l ikely to arise near the shell 
of a hot (10 6 K) highly ionized, w ind-driven bubble. lYorke et al.l (119841 ) used the champagne 
flow model (iTenorio-Tagld Il979l ) - disruption of neutral clouds adjacent to Hll regions- to 
explain the broad lines. The gravitational model (jTerlevich fc Melnicklll98ll ) is based on the 
assumption that the dominant influence is the gravitational potential. In this model, the 
observed Gaussian line profiles are a direct consequence of the virialized motions of the gas 
plus stars complex. If one supposes a self gravitating system, the standard deviation of this 
Gaussian profile in velocity is 



a 



GM 
~R~' 



The Cometary Stirring Model (CSM) proposed later by iTenorio-Tagle et al.l (119931 ) im- 
proves the simpler original picture. They propose that a group of pre-main-sequence low- 
mass stars of a recently formed starburst, move under the gravitational potential with veloc- 
ity dispersion a star - These stars drive winds that stir the remaining cloud through supersonic 
bow, or "cometary", shocks, providing it with an average turbulent motion a gas ~ a sta r- This 
model estimates the total power of the shocks and associated line the lumi nosity which agree 



with the observed rel ation L ~ a 4 . Few attempts have been made (e.g. lOstlin et al.l 12004 



Cumming et al.ll2008l ). and none with a statistically significant sample, to definitely test the 
hypothesis of a gas = cr star which would strengthen the support for this model. 



A skeptical view of this scenario was presented by IChu fc Kennicuttl (119941 ) in a de- 
tailed study of the GHnR 30 Doradus in the Large Magellanic Cloud (LMC). They showed 
that a smooth integrated velocity profile in the inner region (135 x 135 pc) results from 
a complex velocity field, which demonstrates the futility of using global velocity profiles to 
infer the physical origin of the gas motion. The dominant contributor to the global velocity 
dispersion, as proposed by the authors, at least in the bright core, is the superposition of 
individual expanding shells. Consequently, the most important physical mechanisms would 
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be the stellar winds from OB associations combined with champagne-type flows and possibly 
supernovae. 

This long disagreement in the literature has also been enlightened by works favored 
by hig h spectral, spatial as well a s bi- dimensional spectroscopy of GHllRs and the nearest 
HuGs. iMunoz-Tunon et all Jl996h analyzed the velocity field of NGC 604 and NGC 588 in 
the galaxy M33. They show that, in both cases, the peaks of high velocity dispersions are 
systematically concentrated in regions of faint nebular emission. On the other hand, there is 
a smooth low velocity dispersion component which permeates all regions of the nebula. They 
proposed a simple criterion to identify expanding shells in the intensity vs. a diagram. This 
smooth low velocity dispersio n component is what one measures in the integrated spectrum 
of a GHllR. lYang et al.l (119961 ) came up with basi c ally t he same conclusions in a similar study 
of NGC 604. The results found by iTelles et al.l (120011 ) using Echelle long-slit spectroscopy 
of a small sample of HuGs point in the same direction. Despite the complex velocity fields 
found in such objects, they show that the brightest knots dominate the global luminosity 
and velocity dispersion. 

Though much effort has been spent to understand the kinematics of the gas-stars com- 
plex in GHllRs and HuGs, a clear picture is still absent, given the complexity of these 
systems. The puzzle can be noted in the follow question: Why should one find such para- 
metric relations in these violent environments where massive star evolution seems to play a 
significant, maybe dominant, role? Here, we aim to bring some insight to this question by 
extending the bi-dimensional spectroscopy study of Local GHllRs to this nearby prototypical 
HnG, with comparable spatial and spectral resolution. The kinematics of the ionized gas in 
the very center of II Zw 40 is investigated in some detail. 

II Z w 40 has coordinates R .A. 05 h 55 m 42 s . 6, DEC. +03° 23 m 32" (J2000.0), a V = 
15 mag (ITelles fc Terlevichl 119971 ). It lies at a low galactic latitude of ~ —10 degrees. The 
observed radial velocity of 778 kins" 1 and a cosmological (corrected from heliocentric to 
cosmic microwave background frame) recession velocity of 845 kins" 1 , puts it at a distance 
of 11.9 Mpc0. 

The paper is organized as follows. Section [2] describes the observations and reduction 
procedures. The results are presented in § [3] together with full descriptive data analysis. In 
§ H] we discuss our results. Finally, we give a brief summary in § [51 



^Hq = 71 kms 1 Mpc 1 is throughout this paper 
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2. Observations and Data Reduction 

The observations of II Zw 40 were done during the run GN-2003B-Q-26 on 2003 Decem- 
ber 26. We used GMOS equipped with an Integral Field Unit (IFU) on the Gemini-North 
Telescope. The IFU was used in one-slit mode giving a 3.5 x 5 arcsec 2 field of view. Each 
lens in its original hexagonal shape has 0.2" nominal size projected in the sky. In this case, 
the IFU directed to the science target has 20 x 25 lenses, in a total of 500 lenses. A set of 
249 lenses, 1 arcmin away from the science target, provides a sky background for subtraction 
during the data reduction. In order to cover a larger field, we exposed at six different offsets 
overlapping the fields about 1". Figure [T] shows the ACS/HST high resolution Ha image 
(0.027 arcsec/pixel), with the different fields observed with GMOS/IFU shown as rectan- 
gles over the image. The GMOS/IFU data were recorded in three 2048 x 4608 CCDs with 
13.5/im pixel size. We used the R831 grating covering the wavelength range 5970-8080A 
and providing a spectral resolution of 1.5 A (FWHM) at 6580 A, with corresponding Oi ns t = 
29. Oil. km s _1 (er=FWHM/2.355). The exposure time of each field was 300 seconds. The 
seeing remained constantly sub-arcsecond during the observations (~ 0.5-0.6"). 

The basic data reduction was performed using the Gemini IRAfI package developed 
by the Gemini Observatory staff. A CuAr spectrum was obtained in order to provide the 
wavelength calibration. We processed bias subtraction, dome and twilight flat field correc- 
tion. Each field was calibrated in wavelength but not in flux. The final cube for each field 
was re-sampled to 0.2 arcsec per pixel using GFCUBE, resulting in a 16 x 25 x 6270 data 
cube. We developed a special procedure using IRAF in order to produce the Ha emission, 
radial velocity and velocity dispersion synthetic maps. For each data cube we used the task 
FITPROFS in onedspec package to fit a single Gaussian profile to the Ha emission line. The 
task provides the Ha central wavelength, the FWHM and the integrated intensity in relative 
counts. We run the task automatically for all spectra in the original sampling (0.2" hexago- 
nal spaxel) and re-sampled (0.2"/pixel) data cubes. There is no difference using either data 
sets, nevertheless, it is more practical for the analysis to use the re-sampled square pixel data 
cubes. We constructed FITS images using Ha line intensities, centers and FWHM values to 
build our final maps. Finally, we overlapped them to build the mosaiced field. 

As an additional test, th e different maps were also produced from the original data cube 



using the ADHOC package (jBoulesteixi Il999l ) in a parallel reduction procedure in order to 



check the accuracy of the data reduction. This package is designed to reduce Fabry-Perot 



*IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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interferometry data and is able to handle any data cube. Both reductions showed very 
similar maps, then we decided to use the maps produced with the IRAF routines for all our 
further analysis. 

Figure [2] shows a spectrum of a single pixel on the brightest region on the main knot of 
II Zw 40 nebula (knot A). This one pixel spectrum illustrates the quality of our data, and 
shows that we are able to detect, resolve and measure faint lines of [OI]A6300, [SIII] A6312, 
[NII]AA6548,6584, HeIA6678, [SII]AA6717,6731, [NeI]A7024, HeIA7065, [ArIII]A7136, HeIA7281, 
[OII]AA7318, [OII]AA7330, and [ArIII]A7751(not shown). The peak intensity of Ha in this 
one single pixel exceeds 7000 counts and with a total flux of 20000 counts. The S/N on the 
continuum around Ha (which is not of interest here) is about 5 on this one pixel spectrum. 
On the fainter regions, where the line widths are still measurable, the S/N on the adjacent 
continuum is about 0.5, and what determines the accuracy of the line width measurement 
is the peak intensity over the continuum which is directly proportional to the total flux of 
the emission line. We estimated, therefore, a conservative minimum relative flux necessary 
to accurately measure a single Ha line width of 300 counts (log F = 2.5). 

The single Gaussian fit is an operational measurement procedure. Though the major 
contribution comes from regions presenting line profiles well fitted by single Gaussians, some 
regions of the nebula show asymmetric line profiles or multiple components. These measure- 
ments can map regions of kinematic interest through line width variations, and can reveal 
large scale motions through variations of the line centroids across the nebula. When one has 
sufficient high spectral resolution the analysis of a single line profile can provide additional 
information constraining the larger and the smaller scale velocity variations. 



3. Results 

3.1. Electron density and abundance distribution 

We briefly comment on the distribution of some derived physical properties before we 
proceed with the main topic of the kinematic properties of the starburst region in II Zw 40. 
We have been able to reliably measure and map the distribution of some weak lines in the 
central field. 

Figure [3] shows the electron d ensity (n e ) distribution derived from the ratio of [SII]AA6717,6731 



lines ( iDopita &: Sutherland! 120031 ) . This plot was produced by azimuthally averaging the 2D 
map centered on the peak distribution of the sulfur line ratio which is coincident with the 
peak Ha line emission from knot A. Firstly, we show that n e is not constant and ranges from 
~1400 cm -3 in the center to ~ 500 cm -3 about 100 pc from the center of knot A, and higher 
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than the low density limit of 100 cm -3 . It is also noticeable the good fit to a r 1 / 4 law (solid 
line) to the outer points (avoiding the dominant effect of seeing in the center represented 
by the dotted line), which is usually representative of the surface brightness distribution of 
stellar systems in dynamical equilibrium. 

Figure H] shows the oxygen abundance distribution in this central field. Ha and contin- 
uum emission are also plotted as contours over the image. All pixels with Ha emission less 
than 3.4% of the peak e mission were erased. The O/H abundance was derived by applying 



the N2 calibrator from iDenicolo et al.l (120021 ) . This is an empirical relation between the 



the line ratio of [NII]A6584/Ha with the oxygen abundance. Simple closed box chemical 
evolution models considering the primary and secondary origin of nitrogen predict this lin- 
ear correlation between nitrogen and oxygen, though dependent on the ionization parameter. 
Here, we simply apply the empirical calibration of this relation. We found a lower abundance 
over knot A, 12+log(0/H) = 7.85, a slight increase over knot B, 12+log(0/H) = 8.0 and 
a significant increase in a region at NE of the knot A, 12+log(0/H) = 8.3. The statistical 
error of the N2 estimator is $(o/n) ~ 0.15, however relative differences between the regions 
should be real. 

We will not explore the physical conditions analysis here any further. However, we may 
say that we have not found any relation between the velocity dispersion and the oxygen 
abundance or n e . 



3.2. The Ha Integrated Line Profile 

We now compare the integrated line width for the main starburst region in II Zw 40 
derived here, w ith the inner core l ine width obtained with the high dispersion fiber spectro- 
graph FEROS fjKaufer et al.lll999l ). on the 1.52m telescope (ESO]j. We combined 181 specti 



era 

from our GMOS data cube, around the inner core in a synthetic aperture of a diameter of 
2.7" (equivalent to 150 pc) which corresponds to the FEROS fiber diameter. 

Figure [5] shows the line profiles, their single Gaussian fits and the velocity dispersions 
derived from the data of these two very different instruments. Our simulated spectrum from 
GMOS/IFU data shows an observed FWHM = 2.34A, while the FEROS spectrum shows 
FWHM = 1.83A. After the respective instrumental broadening, and the thermal broadening 
corrections (for T e = 10 4 K, FWHM tt = 0.47A or a th = 9.1 kins" 1 ), the velocity dispersions 



§ All the FEROS data and analysis were presented in V. Bordalo's Ms. Sc. dissertation (2004) directed 
by Dr. E. Tellcs at Observatorio Nacional - MCT/Brazil, and will be published in a forthcoming paper. 
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derived were virtually the same, o = 34 kms 1 . 

II Zw 40 presents an integrated supersonic line width with a derived high velocity 
dispersion reaching 34 kms -1 , which is signific antly higher than 30 Dor with a =26 kms " 1 



the most powerful GHllR of our Local Group, (IChu fc Kennicuttlll994l ; iMelnick et al.lll999l ). 



3.3. The Ha Flux Monochromatic Map 

The Ha flux monochromatic map is a result of the emission line intensity measurements 
across the observed field. We fitted single Gaussians to the Ha line profiles to derive their 
respective total fluxes, using the procedure discussed in Section [2j The Ha monochromatic 
(total flux from the single Gaussian fit) contour map is presented superposed on the velocity 
dispersion and radial velocity maps (Figures [6] and [7]) in six contours levels (68, 20, 9, 5, 3.4 
and 2.6% of the peak intensity). The overall Ha distribution in the synthetic map obtained 
from GMOS/IFU agrees very well with the narrow band Ha image from ACS/HST at higher 
resolution (Figured]). In the HST's image it is possible to clearly resolve the two inner knots 
(A and B), the cavities formed by low Ha emission maybe associated with superbubbles, 
and also some filaments. A synthetic red continuum image from our IFU data, shown here 
only in two contours in Figure HI also reveals the two inner knots, separated by projected 
distance of 1.6" (~ 90 pc). However, the Nort hern knot (A) is r esponsible for virtually all 



the ionizing luminosity in the core of II Zw 40 (IVanzi et al.ll2008l ). 



The starburst region revealed in II Zw 40, from the morphological point of view, is also 
very similar to the largest and extensively studied counterparts in the Local Group, such as 
30 Dor, NGC 588 and NGC 604, showing complex kinematic features, along with filamen 



tary structure, shells, and cavities, suggesting a common physical cause (jTenorio-Tagle et al. 



20061 ). The similarities of the properties of the ISM in GHnRs and Hl lGs were first inves- 



tigat ed through empirical relations of integrated physical parameters (IMelnick et al.lll987 



19881 ). particularly the fact that both exhibit supersonic motions of the warm interstellar 
medium. A more detailed analysis of spatially resolved properties of these structures will 
help us understand the interplay between the massive cluster formation and evolution with 
its surrounding interstellar medium on scales of tens to hundreds of parsec. 



3.4. The Ha Velocity Dispersion and Radial Velocity Maps 



The velocity dispersion map is presented in Figure El with Ha intensity contours su- 
perposed. The region corresponding to the brightest knot (inner 50 pc) presents a low to 
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intermediate a value (25-35 kms -1 ). There are, at least, three well defined high dispersion 
(45-65 km s _1 ) regions ~100 pc from the brightest knot. An almost constant velocity disper- 
sion field (< 35 kms" 1 ), which contains most Ha contour levels, is present in a region with 
100x100 pc 2 at NW of the brightest knot. At least two regions of low velocity dispersion 
(15-25 kms -1 ) are also found in the field, strikingly in between regions of high dispersion. 
A third low dispersion region is seen at extreme West, 300 pc far from the core, but we 
preferred to discard this from the analysis due to its low S/N. 

The observed Ha radial velocity map is shown in Figure [71 There is no sign of an 
overall rotation in the main starburst region of II Zw 40 ( 500x500 pc 2 ). On the other hand, 
discrete regions can reach high values of radial velocity, suggesting local expanding motions, 
although the total velocity range does not exceed 60 kms -1 . 



3.5. Intensity-a and Intensity-V elocity Diagrams 



An efficient method to identify kinematic features in GHnRs, detect the massive clusters, 
their effect on the ISM, and their star formatio n rates, is by the analysis of the Intensity-a 



(I -a) and Intensity-Velocity (I-V) diagrams (IMunoz-Tunon et al.l Il996l ; lYang et al.1 Il996 
Martmez-Dekado et aDl2007h . 



We present, in Appendix |Aj a schematic illustr ation of these diagrams, an d a brief 
description of how they were previously interpreted by IMunoz-Tunon et al.l (119961 ). We also 
extend the analysis of these diagrams, including V-a, to bring additional insight to the 
observed bulk motions. 

In Figure [8] (left panel) we show the I-a plot for the whole observed field. The plot 
for II Zw 40 shows exactly the same basi c features found in previous work on GHnRs 
(IMunoz-Tunon et al.lll996l ; lYang et al.lll996l ). In general low intensity regions have the high- 
est a values while the range of a values decreases for high intensity regions. In the case of 
II Zw 40, there is a horizontal band defining a lower limit of cq — 25 km s _1 for regions 
brighter than 4000 counts (log F = 3.6). The horizontal line in the I-a plot shows the veloc- 
ity dispersion of 33.5 km s _1 derived from the integrated line profile with the 2.7" aperture 
(Figure [5] left panel). It is a char acteristic value for I I Zw 4 found in all intensity regions. 
The inclined bands identified by IMunoz-Tunon et al.l (119961 ) are also identified in the plot. 
High velocity dispersions reach peaks of 65 km s _1 in low intensity regions. As we will show 
in the next section, two inclined bands in this plot are superposed in the range 2.6 < log F 
< 3.1. 



Some characteristic line profiles are also shown in Figure [H] (left panel). These are 
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individual pixel profiles exemplifying the variety of shapes observed in different regions. 
Profiles (a) and (b) are the most asymmetric profiles observed. They originate close to 
the center of the two regions of high dispersion at SE and at SW of the brightest knot, 
respectively (see Figure E]). They are clearly double-peak profiles, consistent with expanding 
motions. Profile (c) originates in the second brightest knot and shows a prominent red wing. 
Profile (d) originates in the brightest knot and profile (e) comes from the lowest dispersion 
region. 

Figure M (center panel) presents the I-V plot which shows the relative radial velocity 
variations. These motions are usually associated with the line broadening component named 
(J ra( i- The physical mechanisms at play which produce these line centroid variations should be 
several, such as turbulence, winds, champagne flows, expanding bubbles or rotation, acting 
on stellar radius scales to hundreds of parsecs. However, the estimate of the contribution of 
these subcomponents is limited to the spatial and espectral resolution of the observations. 
For instance, in this work we chose to associate the action of the spatially resolved bubbles, 
producing the highest radial velocity variations, to the component a ra d- The other mentioned 
mechanisms, however, should also contribute to produce the overall pattern in the I-V plot. 
Another line broadening component named u s heii will be here associated with the unresolved 
expanding shells identified as responsible for an underlying broad component found in the 
most inner region and discussed in section 13.71 Therefore, the position-to-position motions 
found in Figure [HI (center panel) contribute to the total broadening of an integrated line 
profile of the whole region. We do not intend to quantify precisely the contributions of the 
broadening components from the diagnostic diagrams. In the case of II Zw 40, the radial 
velocity differences are not higher than 15 km s -1 in the brightest central region, and the 
highest differences are concentrated in low Ha intensity regions. We found that the discrete 
tails toward high intensity values shown in Figure [8] (center panel) are associated with the 
presence of knot A (main tail) and B with a relative radial velocity of ~ -13 km s -1 . 

Figure [8] (right panel) shows the V-a plot. In the case of II Zw 40, it suggests that 
the overall motions in the star forming region are primarily random, instead of dominated 
by systematic line of sight motions. For instance, if one was observing a highly collimated 
outflow toward the line of sight, as a champagne flow, one should expect to see a trend in 
this diagram. Broad profiles from the low density regions should be found blueshifted, while 
narrow profiles from high density regions should be found redshifted, producing a pattern 
such as a narrow band with negative slope, not found for the whole region in Figure [8] (right 
panel). On other hand, the most important result shown by this diagram is that low a (15-25 
kms -1 ) regions have small radial velocity variations, suggesting that they are close to the 
rest frame of the galaxy. This is the first indication that these regions may be unperturbed 
by massive stellar evolution. 
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We will try to associate the morphological features observed in the Ha flux , velocity 
and dispersion maps with the features from these diagnostic diagrams by choosing individual 
regions for a more detailed analysis below. 

3.6. The Kinematics of Individual Regions Revealed by Diagnostic Diagrams 

The patterns shown in I-a, I-V and V-a plots (Figure [8]) re-display many of the features 
that can be observed in the maps (Figure [6] and Figure [7]). We wish to evaluate the contribu- 
tion of individual features to the integrated Ha line profile of II Zw 40 using the kinematic 
diagnostic diagrams. For this purpose, we defined eight regions with reliably measured a and 
radial velocity. The choice of regions was arbitrary, but motivated by different peculiarities 
noted in the Ha flux, radial velocity and velocity dispersion maps. The chosen regions are 
superposed on these maps in Figure [91 

For each chosen region, we have plotted a set of four graphs. These graphs are presented 
in Figure [10] in blocks of four boxes: (a) an integrated profile over the region in linear- 
logarithmic axes, where we also show the region number, the derived a value, and the 
estimated skewness £ and kurtosis n of each integrated profile, (b) I-a, (c) I-V , and (d) the 
V-a plot. The vertical solid line in plots (b) and (c) shows our confidence limit in flux as 
described in Section [21 All references to the boxes in this section refer to Figure [101 where 
we quote the region number and the corresponding box letter, i.e. box (la) means plot (a) 
of Region 1 in Figure [10J The results of the analysis of the emission line shapes through £ 
and k estimates for each region will be presented in Section 13.71 A brief description of each 
region follows: 

Region 1 was defined to cover the main and brightest Ha knot (knot A) as shown by 
the Ha map (Figure [9] left). It covers an area of 1.8" x 1.6" (equivalent to 100 pc x 90 
pc). The resultant integrated profile of Region 1 is shown in box (la), where the dashed 
and dotted lines represent the instrumental profile and the Gaussian fit, respectively. The 
a derived for Region 1 was 33.3 km s -1 which agrees with the value derived from our IFU 
simulated and FEROS spectra for a larger aperture (see Figure [5]). Box (lb) shows that there 
is a well defined lower cq value slightly above 25 km s -1 and no high a variations. Radial 
velocity variations found in box (lc) indicates that the gas in Region 1 expands at < 15 km 
s _1 . Since the points in box (lc) are within ~ 7.5 km s _1 from the average recession velocity, 
it is likely that the rms radial velocity variation contributes to an equivalent FWHM of ~ 15 
km s _1 or a ~ 6.5 km s _1 (see Appendix [Aj. Radial velocity variations between individual 
pixels in Region 1, therefore, do not contribute significantly to explain the total 33 km s _1 
velocity dispersion. At this spatial resolution only large scale radial velocity variations can 
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be identified, and only these will be analyzed here (~ 50 pc). Box (Id) shows that most 
points are clustered with no apparent correlation between V and a. We found that random 
motions are dominant in Region 1, instead of line of sight motions. 

Region 2 was defined to cover the second brightest knot (knot B), south of Region 1. 
Region 2 covers an area equivalent to 90 pc x 90 pc. Box (2a) shows that, the integrated 
profile is the most asymmetric of all regions, showing a prominent red wing. This region is 
responsible for the red wing in the integrated profile obtained for II Zw 40 with FEROS, and 
also with our GMOS/IFU data (Figure The reason for this asymmetry is that a reaches 
a local maximum in a few relative high intensity points (North of knot B) as shown by box 
(2b). Moreover, the same high intensity points, seen in box (2c), have the largest recession 
velocities, resulting in this observed red wing. As a consequence, box (2d) shows a strong 
correlation pattern (see Figure [T3"C). The radial velocity motions in Region 2 contribute 
twice as much to the integrate line width of the region as compared with Region 1 (a ~ 13 
kms" 1 ). The gas in Region 2 has a relative line of sight motion with respect to the gas in 
Region 1 of roughly -13 kms" 1 . We found, however, two systematic components, one due 
to the relative radial motion between the two knots, and the other due to the expanding 
motion, as shown by the vertical band at log F ~ 3.2 in box (2c) presenting an intermediate 
radial velocity variation. 

Region 3 covers the SW high dispersion region. It covers an area of 1.8" x 1.6" 
(equivalent to 100 pc x 90 pc). The expected kinematic features due to the presence of 
an expanding shell are seen in all boxes for this region. Box (3a) shows a symmetric broad 
profile with a = 45.4 kms -1 . This is actually the resulting smoothed profile obtained due 
to the contribution of unresolved asymmetries and multiple components as shown in one 
pixel profiles (a) and (b) of Figure [8] (left panel). Box (3b) clearly shows the inclined band 
indicating the presence of an expanding shell, possibly produced by an expanding bubble (see 
Figure [TBIA). Box (3c) confirms the expanding motion through the very high spread in radial 
velocity (see Figure [T3B). The observed range of radial velocity, of about 20 kms -1 above 
and below the average recession velocity should underestimate the real radial component. 
This is in fact a limitation of fitting single Gaussians to the integrated profiles. From the 
line splitting in individual pixels in the center of the bubble we could measure a maximum 
AV ~ 90 km s -1 or a V exp ~ 45 km s _1 . We estimated the radial velocity flow component 
to the integrated profile shown in box (3a) as FWHM ~ 90 km s" 1 or o ~ 38 km s -1 (see 
Appendix [A]). Therefore, the expanding motion dominates the total broadening seen in the 
integrated profile of Region 3. The data points in Box (3d) show large scatter and there is 
not a significant trend in this plot, as seen in Region 2. 

Region 4 covers an area of 1.8" x 1.8" (equivalent to 100 pc x 100 pc) and presents 
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a field of intermediate a and Ha intensity values, north of knot A. Region 4 presents many 
similarities with Region 1: similar a derived as shown in box (4a) and similar pattern 
in box (4b) showing low a variations and a lower limit ~ 25 kms -1 . The small radial 
velocity variations found in box (4c) and a concentrated velocity dispersion distribution in 
(4d) indicate that there is no significant systematic motions contributing to the total line 
broadening, which is dominated by random motions on scales smaller than the seeing of ~20 
pc. 

Region 5 is one of the most interesting regions, and covers a large area of 2.6" x 2.6" 
(equivalent to 150 pc x 150 pc) SE of knot A. It presents some of the broadest profiles, and 
it is probably associated with an expanding shell (boxes 5a and 5b). Region 5 also presents 
high radial velocity spread (box 5c) of over 40 kms -1 , although this may underestimate the 
real radial velocity component as argued above in the analysis of Region 3. Regions 3 and 5 
show very similar patterns in all diagrams, suggeting that these patterns can be easily used 
as expanding shell diagnostics in GHllRs. Despite their different sizes and morphologies 
the two shell regions have similar kinematic properties. Some double peak profiles seen in 
Region 5 are also consistent with a bubble expanding with a velocity of 45 kms -1 , likely to 
be powered by supernovae and stellar driven winds. There is not a significant trend found 
in box (5d). 

Region 6 was defined to cover a field containing the lowest velocity dispersion values. 
It covers an area of 2.0" x 2.0" (equivalent to 110 pc x 110 pc) North-East of Region 1. The 
observed a = 23.5 kms -1 , as shown in box (6a), is the lowest value verified of all regions. 
Box (6b) also shows that a significant area presents a values lower than knot A, though they 
are still supersonic. Velocity variations are also very small as shown by box (6c). Box (6d) 
indicate primarily random motions. All plots suggest, therefore, that Region 6 represents 
a region where the mechanical energy input resulting from massive star evolution does not 
play a significant role to stir the gas up. 

Region 7 presents a relatively low velocity dispersion field (box 7a). It covers an area 
of 1.6" x 2.0" (equivalent to 90 pc x 110 pc), South of knot B. It completes the sampling 
of a larger low-to-intermediate dispersion field covered partially by Region 2 (boxes 2b and 
7b). Both regions are confined between the two high dispersion bubbles, Regions 3 and 5. 
The Ha structure in Region 7 is associated with the extended emission to the South of both 
knots (see Figure [9] left panel). Box (7c) shows a weak correlation of / and V with fainter 
emission points moving away. However, no clear correlation between radial velocity and a is 
seen in box (7d). 

Region 8 covers an area from W to NW of the knot A (equivalent to 70 pc x 160 
pc) with a cr = 32 km s -1 (box 8a), similar to Region 1. It was defined to contain a field 
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with almost constant velocity dispersion through a large range of intensities, which is shown 
by box (8b). The pattern formed in box (8b) is not an inclined band. Region 8 does not 
present high velocity variations, as shown by box (8c) and it is likely dominated by random 
motion inferred by the points concentrated in box (8d). The few higher dispersion points at 
low intensities seen in box (8b) and (8d) may be picking up a possible additional region of 
interest (i.e. shell) North of Region 8 (see Figure [9] right panel) which is not well mapped in 
our observations and is not included in the present analysis. 

Our definitions of the regions favor the interpretation of the different patterns in I-a 
diagrams where kinematic features are easily identified. Thus, the patterns may be the 
result of different physical mechanisms. Something to have in mind here, is our spatial 
resolution which only permits us to identify kinematic features in scales of a tens of parsec 
and limited by our field size. The overall picture is that local systematic motions represent 
a dominant component to the observed line widths in regions 3 and 5 which are likely 
associated with resolved expanding shells due to the fact that we can observe line splitting 
directly. Although, masked by the use of single Gaussians, they leave a signature in the 
diagnostic diagrams. Region 1 presents a small range of intermediate a values which are also 
found in all the intensities of nebular emission across the whole extent of the star forming 
region. Region 2 is turbulent with some characteristics of flow motions. Regions 4, 7 and 8 
show smooth turbulent fields. Wherever we cannot observe double peaked profiles in these 
regions, we cannot completely rule out a radial systematic component to the line width, 
possibly associated with unresolved shells, however random motions must dominate over 
expansion in these cases. Finally, relatively to the other regions, Region 6 seems to be 
slightly affected or even unperturbed by mechanical energy input provided by massive star 
evolution, presenting the lowest a values. 



3.7. Making up the Integrated Velocity Profile from Individual Regions 

Figure [TT] shows the fully integrated line profile produced by combining all eight spectra 
for the regions presented in § 13.61 with a resulting a = 34.1 kms" 1 . Despite the complex 
medium-scale structures (50-100 pc) seen in the synthetic maps (Figures [6] and [7]), mainly 
due to the presence of asymmetries and multiple components (e.g., profiles (a) and (b) of 
Figure [HI left panel), the fully integrated profile is quite smooth and symmetric. 

An important result is that the velocity dispersion derived from a single Gaussian fit 
to the profile of Region 1 (Figure [10] box la) is virtually the same as those derived from the 
original FEROS spectrum with an aperture of 2.7", and the simulated spectrum with the 
FEROS aperture on our GMOS/IFU data (Figure ED, and the fully integrated line profile 
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(Figure [TT]). Region 1 represents 10% of the total area analyzed but 49% of the total Ha 
flux in all eight regions. In contrast, Region 3 and 5, showing the broadest profiles, cover 
33% of the total area analyzed but only contribute 14% to the total flux. Other regions, 
namely Region 6 and 7, covering 25% of the total area and showing the lowest o values only 
contributes 10% to the total Ha flux. All the highest surface brightness regions, namely 
Region 1, 2, 4 and 8, show the common values a = 32-35 km s _1 which characterize the 
global kinematics of the warm gas in II Zw 40's dominant starburst region. 

The aperture effect to derive the characteristic velocity dispersion of II Zw 40 is negligible 
if an int egrated spectrum is taken with an aperture covering the brightest knot. As an 



example, iMelnick et al.l (119881 ) derived for II Zw 40 a a = 35.2±0.5 km s _1 with a 6" wide 
entrance slit aperture. From an observational point of view, this result is very positive. If the 
same is valid for most HnGs and GHnRs, as has been confirmed, the measured supersonic 
a, introduced in the L-a relation, is little affected by the effects of SNRs and wind driven 
bubbles where the warm gas emission is faint, and where the integrated velocity dispersion 
is much higher than in the central core. 

We computed the skewness, £, and the kurtosis, or flatness factor, k, to quantify the 
shape of each profile and they are included in all boxes (a) of Figure [10J We warn that 
the quantitative results, especially for k, are sensitive to the window of summation and the 
height of continuum. We defined a confident fixed window 6576-6584A to compute the £ and 
k values, since broad wings are not well mapped for the faintest profiles. Relative differences 
between kurtosis values are more significant than their absolute values, which should in 
principle be comparable with a Gaussian profile (£=0, k=3). We found that most profiles 
have peaked shape k > 3 (leptokurtic). Region 3 and 5 are significantly more flattened 
than the others (platykurtic). The most asymmetric profiles identified by the highest | £ | 
are those of Region 2 and 6, which can be easily distinguished visually in liner-logarithmic 
plots in boxes (2a) and (6a) of Figure [TUJ We also estimated £ and k values for the fully 
integrated profile shown in Figure [TT1 We can see that not only the a but also £ and k values 
are virtually identical to those estimated from profile of Region 1 (Figure [TOl box la). 

Since the S/N in the wings of the profile from Region 1 is much higher than the others, 
we were able to measure across in a larger window (6574-6586A) and estimate another pair 
of £ and k values for this profile, giving an idea of how sensitive the values of the higher 
order moments are. The estimated values derived in this spectral range are £=0.14 and 
k=4.19, showing that the integrated profile in Region 1 is very peaked. This result indicates 
the presence of a broad but weaker component in the core of II Zw 40. Multicomponent 
Gaussian fits to the profile of R egion 1 are shown in Figure [121 We used the PAN routine 
(Peak ANalysis; Dimeo ( 20021 )) in IDL to fit two components to Ha profile. We found a 
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narrow component a n = 29.4±0.6 km s _1 which is not very different from a derived from a 
single fit a — 33.3±0.9 km s _1 . The broad component was estimated to be cr& = 73.2±1.0 
km s _1 and its flux corresponds 18% to the flux of the narrow component. We may speculate 
that the broad component is associated with unresolved wind-driven shells at the core of II 
Zw 40 where most of the young and massive stars concentrate. 



4. Discussion 



An important verification provided by the above results is that the Ha kinematic features 
already found in Local Group's GHnRs are represented in the main starburst region of II 
Zw 40: a prototypical Hn galaxy at roughly 11.9 Mpc. In fact, one way to investigate the 
starburst phenomena on these galactic scales is to compare the main starburst region in II 
Zw 40 with the local GHnRs studied by similar techniques. 



Yang et al.l (119961 ) found that virial motions and expanding shells contribute roughly 



equally to the veloci t y widt h of the integrated profile in NGC 604. Using the same methods, 
Muhoz-Tuhon et al.l (119961 ) proposed an evolutionary scenario when they compared their 
results of NGC 604 with NGC 588. A simple model for an idealized homogeneous shell 
evolving in a uniform density medium suggests that younger and fast shells should have 
lower intensities and reach higher a values than older ones. Hence, the authors argue that 
NGC 604 appears to be younger than NGC 588 based on the inclined bands associated 
with shells, and identified in l-a plots. Furthermore, NGC 604 presents a well defined 
lower limit (o"o — 17kms _1 ) in supersonic line width everywhere (horizontal band in l-o 
plot), while NGC 588 does not. A substantial fraction of NGC 588 presen ts line widths 



l ower than the brightest region, including points showing subsonic a values. iMelnick et al. 



( 119991 ) suggest that it is also the case for 30 Dor. In the Cometary Stirring Model (CSM) 
framework, this difference can also be understood as an evolutionary effect. A young region 
keeps most of the supersonic motions communicated to the leftover gas via bow shocks 
by the winds of low-mass stars during its early stage. As massive stars evolve, ejecting 
mechanical energy into the ISM, through stellar winds and supernova events, they produce 
shells that expand beyond the core. Shells will be destroyed by random motions of the 
background neighboring gas and slow down, resultin g in a even below the bright region 
( jTenorio-Tagle et al.l Il993l ; iMunoz-Tunon et al.lll996l ). The main point in this scenario is 
that the dominant line-broadening mechanism in young GHnRs is gravity, while tuburlence 
by wind-driven shells and supernovae dominate in older ones. 



More recently, iMartmez-Delgado et al.l (120071 ) have shown that the morphological fea- 
tures observed in their analysis of three-dimensional spectroscopic data of three blue compact 
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galaxies (Mrk 324, Mrk 370, and III Zw 102), especially those in the a versus I pea k diagram, 
look similar to the ones found in GHnRs, with a supersonic horizontal band that extends 
over a large range of intensities, and inclined bands that reach high a values at low intensity 
arising from large volumes surrounding the centr al knots. In a simila r study to our own, 
using similar technique with Gemini GMOS/IFU, Westmoquette et al. ( 2007al lbh looked at 
the properties of a young star cluster and its environment in the dwarf irregular starburst 
galaxy NGC 1569. Despite the fact that they decompose the profiles rather than fit a single 
Gaussian, the narrow components of their profiles show all properties found in these other 
studies including ours, namely that, it is most likely explained by a convolution of the stir- 
ring effects of massive star evolution and gravitational motions. In addition, they found that 
the observed lower limit (o"o — 12 km s^ 1 ) of the narrow component is in agreement with 
^GM/R. 

Here, we wish to test the hypothesis that the dominant broadening mechanism of the 
Ha line in the core of II Zw 40 is gravity. 

As shown in Figure [SJ( left panel), II Zw 40 presents the lowest a values (er ~ 20 km s^ 1 ) in 
regions with faint Ha emission. Here we identify Region 6 as the most significant contributor 
to this low a (see Figure [10] box 6b), which appears to present dominantly random motions, 
through its diagnostic diagrams. In addition, Region 6 does not seem to be affected by the 
neighboring superbubble of Region 5. Therefore, this region may still retain the kinematic 
information of the pre-existing underlying turbulent velocity field. The velocity dispersion 
values derived in Region 6 should not be associated with the internal velocity dispersion of 
individual dense molecular cloud cores or clumps (typically less than a few kms -1 ), where 
individual or groups of massive stars form, but instead, they may be the result of the large 
scale turbulent motions of the parent complex of diffuse molecular clouds, which gave rise to 
the present episode of violent star formation. Some of these clouds will form stars in short 
local dynamical times, while others, though dense, may not form stars. We do not intend to 
discuss any detailed model, but simply try to envisage a viable scenario for our conjecture 
of an overall underlying velocity component. Massive star clusters will form in much denser 
fragments (clumps or cores) of the molecular clouds whose properties are not assessed by 
our observations, and require density tracers (e.g. HCN) with high spatial resolution radio 
observations. Though, these clumps must also have relative velocities stirring up the ISM 
under the influence of the common gravitational potential due to the complex of stars, and 
gas, and possibly dark matter. The 12 CO (J= 1 — ► 0) traces the diffuse molecular cloud 
on hundreds of parsec scales. The measured CO line width of a region encompassing the 
entire ionized region in II Zw 40 is AV = 42 km s _1 or a ~ 18 km s~ l ( jTacconi fc Young 



19871 ). somewhat comparable to the low a in the ionized gas. However, these motions seem 



to be completely detached from the global HI velocity distribution which extends up to 5 



-19- 



arcmin (~ 6 Kpc) in the direction of the southern and northern tails (Ivan Zee et al.lll998 



Brinks fc Kleinlll988l ). where the line width is W50 ~ 120 km s 1 (a ~ 50 km s *) and has 



a strong contribution by the large scale rotation pattern of the galaxy on these scales. 

The large resolved supperbubble in Region 5 is very well delimited in a and intensity 
maps and may not have slowed down sufficiently to be disrupted. Moreover, Region 3 
also presents a well defined bubble in a and intensity maps, which is seen as an inclined 
band in the diagnostic plot. Despite being smaller than the bubble in Region 5, they both 
show similar profile shapes, expansion velocities V exp ~ 45 kms -1 , and integrated velocity 
dispersions a ~ 45 kms -1 , possibly indicating similar evolutionary phase and associated 
with the age of knot B. 

Although II Zw 40 contrasts slightly with a simple scenario for a GH11R, in the sense 
that it may suffer from the composite effects resulting from the evolution of multiple stellar 
clusters, we can say that in the CSM context, the GH11R in II Zw 40 is kinematically young. 
Some superbubbles are evolving within the c ore, though it does not present a constant o or 



only "well-behaved" Gaussian emission lines (jMunoz-Tunon et al.lll995l Il996l ). However, the 



presence of surrounding regions showing lower a than the core in II Zw 40 is a sign of youth. 
There are still unperturbed regions preserving kinematic information of the parent cloud 
under the dominant influence of the galaxy's gravitational potential, ionized by diffuse UV 
radiation. In consonance with our conclusions, many studies about this galaxy indicate the 
presence of a very young starburst in its brightest region, as mentioned in the introduction. 

An impor t ant r esolution bias must be of concern here in order to test our hypothesis. 



Melnick et al.1 (119991 ) show that the supersonic integrated line profile in 30 Dor is actually 
due to a superposition of discrete number of clouds, small shells and filaments which are 
resolved on subparsec scales. In these scales, strong winds of WR, O and B stars might be 
responsible for stirring up the gas around the dominant young cluster, imposing a turbulent 
motion higher than the one predicted only by the gravitational potential. The presence of 
a significant number of the very young stars in the core should be responsible partially for 
the difference between the a found in Region 1 and that one found in Region 6. In fact, 
we found evidence of a broad component cr& ~ 73 km s -1 in Region 1, responsible for the 
large wings in the integrated profile, which may be associated with unresolved inner shells. 
However, its contribution to the integrated profile is very small, since only 18% of the total 
Ha flux in Region 1 comes from this component. 

We can now sum up the main reasoning so far. A common framework assumes that the 
observed velocity dispersion (a b s ) is mainly given by the contributions of gravity (cr grav ), 
thermal broadening (a t h), multiple unresolved expanding shells (cr s heii) and systematic radial 
motions (<J r ad)- Here we consider a ra d as being the sum of several phenomena such as SN, 
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expanding bubbles, champagne flows, relative motions between discrete clusters providing 
large scale shear, turbulence and rotation, which cause position-to-position radial velocity 
variations and, in extreme cases, line splitting. We detected and identified cr s heii on ly i n the 
core, i.e. Region 1, as the weak broad component, but possibly present with different weights 
in the whole nebula. Indeed, we found that most profiles have peaked shape indicating the 
presence of faint wings or broad components, with the exception of those found in regions 
of bubbles, i.e. Region 3 and 5. Once corrected for the instrumental profile, the observed 
velocity dispersion of a GHnR or discrete regions inside them is therefore given by 



(Tabs = \/Vgrav + ^shell + a th + a lad 



(IMelnick et al.lll999l ). We fixed the thermal broadening as Oth = 9.1 km s _1 due to the gas 
electron temperature for hydrogen at 10 4 K, causing no significant systematic error in our 
results. The radial velocity variations on scales larger than 50 pc can be identified here. It 
varies from few km s _1 in regions such Region 1, a ra( i ~ 6.5 km s -1 , to a Ta d ~ 38 km s -1 in 
Regions 3 and 5 as inferred from line splitting in the center of bubbles. Region 6 seems to 
be "free from shells", and it can be a good probe of a grav ~ 20-25 km s _1 , posing a lower 
limit of velocity dispersion in II Zw 40. In fact, these values are not very different from the 
narrow component found in Region 1, a n ~ 29 km s -1 , and the lower limit do ~ 25 km s _1 in 
the I — a diagnostic diagram. If indeed, the a measured in Region 6, is unaffected by stellar 
evolution, and the narrow component of the Gaussian fit to the profile of Region 1 gives the 
same low value of a, and additionally the low limit of overall measured a is probing the same 
physical mechanism, we may suggest a common dominant line broadening cause. Gravity 
provides the continuous underlying source of stirring the ISM over the whole star-forming 
region. However, the characteristic o = 32-35 km s -1 found over the whole starburst region 
of II Zw 40 and derived from a single Gaussian fit should be the result of combined sources, 
of the dominant action of gravity, unresolved multiple expanding shells and systematic radial 
motions. 



5. Summary 

The three-dimensional spectroscopic analysis of II Zw 40, using Gemini-North GMOS/IFU, 
yield insights into several astrophysical questions about the structure and evolution of star- 
burst regions, specially those associated with the birth of super star clusters and their impact 
on the ISM. The kinematics of the warm gas inferred from line profiles of emission lines, such 
as Ha, imposes constraints on the theoretical models of star formation and on feedback of 
stellar evolution on the ISM. We have analyzed the kinematic properties of the ionized gas 
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in II Zw 40 using diagnostic diagrams and individual region description. A summary of the 
main conclusions follows: 

(1) Electron densities in the inner starburst region range from 5 x 10 2 to 1.4 x 10 3 cm -3 , 
somewhat higher than the common low density regime for Hn regions of 10 2 cm" 3 . Its 
spatial distribution is not constant and shows a clear decline (compatible with an r 1 / 4 
law) with a peak coincident with the Ha line emission peak. 

(2) Oxygen abundance as inferred by the N2 calibrator in the inner 100 x 100 pc 2 ranges 
from 12+log(0/H) = 7.8 to 8.4, showing lower values 12+log(0/H) = 7.85 associated 
with the dominant site of star formation, namely knot A. 

(3) The aperture effect to infer the characteristic velocity dispersion a = 34 km s" 1 of II 
Zw 40, derived from a single Gaussian fit, is negligible. All shape properties (a, £ and 
k) of a fully integrated profile are virtually the same as those found in the profile of 
the brightest region. 

(4) The diagnostic diagrams of / versus a, I versus V, and additionally a versus V are 
powerful tools to identify the origin of internal motions in GHnRs and also in HllGs. 

(5) We found regions where stellar evolution does not seem to play a significant role stirring 
the gas up, and should reveal the kinematic signature of the proto-cloud that gave rise 
to the present starburst. In the case of II Zw 40, Region 6 North-East of the main 
knot poses a lower limit of ~ 23 kms -1 to the global velocity dispersion in the galaxy, 
which should provide a good estimate of cr grav . 

(6) Despite the different sizes and morphologies, the two large shells found in Region 3 
and 5 have the same kinematic properties, with very similar patterns in all diagnostic 
diagrams. Both regions have smoothed integrated line profiles with a ~ 45 kms -1 and 
they seem to expand at V exp ~ 45 kms -1 , contributing significantly with a rad ~ 38 
kms" 1 to the total width of their integrated profile. 

(7) The inner 50 pc and brightest region of II Zw 40 seems to be kinematically very young, 
where gravity seems to play an important role. Part of the difference found between 
a ~ 33 km s" 1 in the core and 23 kms -1 found in Region 6 is possibly associated 
with a weak broad component and a small systematic radial velocity flow component 
<y ra d- The broad component should be associated with a s h e ii originated from multiple 
unresolved shells provided by stellar winds of WR and OB stars acting on few or 
subparsec scales. 
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(8) The kinematic features in II Zw 40 are all remarkably similar to the ones found in 
supersonic GHnRs in irregular and other star forming galaxies powered by massive 
star and stellar cluster formation and evolution, suggesting that these starbursts in 
HII galaxies are their scaled up versions on galactic scales. 



Finally, the core supersonic a is the one producing the Luminosity-cr relation observed 
in GHnRs and HllGs. This component can be derived by simply measuring the integrated 
line profiles by fitting a single Gaussian on the brightest knot of the starburst, regardless 
of aperture size. However, we must further investigate how these derived motions precisely 
relate to the underlying mass distribution before we can derive absolute total galactic masses. 
We must also investigate, with a statistically significant sample, the possible evolutionary 
effects of the starburst on the observed relations, and how they can be parametrized for the 
use as a powerful extragalactic distance estimator applied to high redshifts. 
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A. Interpreting the Diagnostic Diagrams 

It is convenient to emphasize the meaning of some patterns formed by the point distri- 
butions in the different diagnostic diagrams of Figure El These patterns are the graphical 
representation of the kinematic signatures observed in the 2D flux, velocity and dispersion 
maps. The following description is valid for the methodology of fitting only single Gaussians 
to the emission line profiles and its limitation is emphasized when necessary. However, it is 
a very powerful tool for kinematic analysis provided by data cubes. Three typical patterns 
are of particular interest, and they are sketched in diagrams (A), (B) and (C) of Figure [T3l 

Figure [T3l\ shows the pattern known as the inclined band in the I -a diagram, previously 



identified by iMunoz-Tunon et al.l (119961 ) . The authors propose that this pattern should be 
the result of an ideal shell evolving in the ISM. Their scenario allows the interpretation of 
this diagram as a diagnosis for the evolutionary state of the GHllRs. The points defining the 
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vertexes of the triangle should correspond to the three parts of the shell: center, inner edge 
and outer edge (see their Fig. 3). The velocity dispersion in the center reaches the maximum 
value, while the minimum value occurs in the outer edge. The inner edge presents higher 
intensity and higher velocity dispersion than the outer edge. Assuming that the pattern 
drawn in the Figure [TBI A represents a young shell, as it ages, a at the center becomes 
gradually lower, as well as the difference between the intensity at the center and at the inner 
edge. Then, the whole figure shifts down, and the upper angle becomes smaller. Furthermore, 
in their scenario, older shells should present higher intensities than young ones, and should 
be mapped at right of the young ones in the diagram. 

Figure [T3B shows a vertical band pattern in I-V diagram. This pattern indicates the 
contribution of resolved systematic radial movements. Abrupt velocity variations in a short 
intensity range should indicate a local expansion of the gas. It is expected that regions which 
present the inclined band pattern in the I-a diagram also present a well defined vertical band 
pattern in I-V diagram. In general, any vertical variations (AV) can be used to quantify the 
component due to expanding motions in a region. We used the approximation FWHM ra( j ~ 
AV to infer the broadening component due to radial velocity variations {cT ra d) for regions with 
almost constant a values and nearly Gaussian emission line profiles. However, if the emission 
profiles are not well represented by single Gaussians, AV measured in I-V plot should 
underestimate the real component due to systematic radial expansion. As an alternative, 
one might try to measure the line splitting from the most asymmetric or double peaked 
profile, in order to better estimate AV from two centroids. 

Figure fTBTJ shows a correlation between V and a. This pattern indicates systematic 
motion with a significant component in the line of sight. The pattern drawn indicates 
differentiated group behavior: in this case, gas with relatively high a is moving away from 
the observer. Although we have drawn the pattern with a positive slope, a region with the 
same systematic motions could be identified in a negative slope pattern. Such dependence 
should be expected in the presence of relative motions between distinct clouds of ionized 



gas with different internal properties. Ch ampagne flows shou 



the one shown in Figure [T3G (see § IIIc in lSkillman fc Balickl (11984 )). The simplest pattern 



d rep roduce the pattern like 



would be the random distribution of the points in the plane V-a, indicating no dependence 
between the variables. In this latter case, it may represent isotropic expansion or turbulence, 
but with no privileged direction. 
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Fig. 1.— A high resolution Ha image (0.027"/pixel) obtained with ACS/HST. The 
GMOS/IFU fields are shown over the central region of the galaxy (HST's arch ive data from 
the proposal 9739 by Rupali Chandar; see also Figure 1 in Vanzi et al. ( 20081 )). 
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Fig. 2. — Single pixel spectrum on the bright core centered on knot A. The labels show most 
of the emission lines present on our spectrum from 6250A to 7400A. The y axis has been cut 
to illustrate the relative intensity of the lines, in particular the much larger intensity of the 
Ha line which is used in this work. 
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Fig. 3. — Electron density distribution in the central field, derived from an azimuthal fit to 
isointensities of the ratio of the sulfur line ([SII]AA 6717,6731) maps. The solid line is a r 1 / 4 
linear fit to the data. The dotted line is a Gaussian profile to represent the approximate 
seeing during this observation. 
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Fig. 4. — Oxy gen abundance distri bution in the central field, derived from the ratio of 
[NII]A6584/Ha flDenicolo et al.ll2002h . The Ha emission contours at 68, 20, 9, 5 and 3.4% of 
the peak intensity (solid lines) and the continuum emission contours (dashed lines) at 47% of 
its peak intensity are superposed on the image. We use a 3x3 box smooth the pixel-to-pixel 
variation due to poisson noise 
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Fig. 5. — Line profiles derived with GMOS/IFU (left) and FEROS (right) data. The dashed 
lines represent the respective instrumental profiles (FWHM = 1.5A for the GMOS and 
FWHM = 0.13A for the FEROS). The dotted lines represent the single Gaussian fits. The 
velocity dispersions values derived are shown in the boxes. 




Fig. 6. — Ha velocity dispersion map, corrected for instrumental and thermal broadening, 
of the inner region of II Zw 40. The contours represent 68, 20, 9, 5, 3.4 and 2.6% of the Ha 
peak intensity. North is up and east is left. 
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Fig. 7. — Ha radial velocity map. The contours represent 68, 20, 9, 5, 3.4 and 2.6% of the 
Ha peak intensity. The recession radial velocity is not corrected by the heliocentric motion. 
North is up and east is left. 
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Fig. 8. — I -a (left), I-V (center) and V-a (right) plots for the whole region observed in II 
Zw 40 as described in the text. The solid horizontal line in I-a plot represents the a value 
for the inner core (D = 2.7"). Characteristic individual pixel line profiles are shown in the 
left panel from (a) to (e). Dashed lines represent the single Gaussian fits and the dotted 
line under profile (d) is the instrumental profile. The vertical solid line at log(Relative Flux) 
= 2.5 represent the estimated minimum line flux to well measure a single line profile. All 
points with log(Relative Flux) < 2.5 are not shown in V-a plot. 
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Fig. 10. — The panel present the four boxes discussed in the text for each region. The 
dashed line in box (la) represents the instrumental profile (FWHM = 1.5 A). The dotted 
lines in all boxes 1 represent the single Gaussian fits to the integrated profiles. Note that 
in linear-logarithmic axes a Gaussian takes the form of a parabola. The region number 
and the a value derived are shown in the boxes (a). I -a, V-a and V-a for each region are 
presented by boxes (b), (c) and (d), respectively. The vertical solid line in boxes (b) and 
(c) at log(Relative Flux) = 2.5 represent the estimated minimum line flux to well measure a 
single line profile. 
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Fig. 10. — Continued. Only Region 6 presents points with log(Relative Flux) < 2.5 and 
they are not shown in box 6b. 
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Fig. 11. — The fully integrated profile summing all eight regions analyzed in section 13.61 
The dotted line represent the single Gaussian fit to the profile. 
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Fig. 12. — The integrated profile of Region 1 (same in Figure fTUI box la) plotted in linear 
axes. The profile is well fitted by two Gaussians and the profiles are plotted. The dashed line 
represents the narrow component a n , whereas the dotted line represents the broad component 
&b- The solid line shows the resultant profile. 
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Fig. 13. — Schematic diagrams showing idealized patterns formed by the distribution of 
points discussed in the text. 



